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Abstract — In this paper, the capacity of wireless channels is 
characterized based on electromagnetic and antenna theory with 
only minimal assumptions. We assume the transmit antenna 
array of arbitrary configuration is confined in a spherical region 
and the receive antenna array is isotropically uniform. The 
capacity is shown as (aP/N ) log e, where aP is the total received 
power and iVo is the noise spectral density. Surprisingly, our 
result shows the capacity is irrespective of the size of the transmit 
antenna array against conventional belief, and is only limited by 
the received power not by the degrees of freedom. 

I. Introduction 

The capacity of multiple-input multiple-output (MIMO) 
channels has been considered to be severely limited by the 
size of antenna arrays. The capacity scaling for a three- 
dimensional network Jl] and the degree-of-freedom analysis 
for a polarimetric antenna array [2] were provided based on 
the spherical vector wave decomposition, and both showed the 
number of usable channels is proportional to the surface area 
enclosing the transmitter. 

There have been a lot of attempts to squeeze more antennas 
into a given space. For example, polarimetric antennas were 
used to get degrees-of -freedom and multiplexing gain by using 
the wave polarization [2], [3|, [4|. Also, the MIMO cube was 
introduced, which consists of twelve dipoles located on the 
edges of the cube to increase the degree of freedom in a limited 
space |5j, (6). Besides, it was recently shown that even when 
the antenna spacing at the receiver is negligibly small, two 
degrees of freedom can be achieved for a two-user multiple 
access channel, which was previously thought impossible (TJ. 

In this paper, we attempt to characterize the ultimate limit of 
wireless communication by deriving the capacity of wireless 
channels with only minimal assumptions. We assume the trans- 
mit antennas are confined inside a sphere of a certain radius but 
otherwise completely arbitrary and assume the total transmit 
power is constrained to be P. We also assume arbitrarily many 
receive antennas so that they do not become a bottleneck. We 
show the capacity is given by (aP/No) log e, where aP is 
the total received power and Nq is the noise power spectral 
density. Interestingly, the capacity is irrespective of the size of 
the source region, which is due the availability of arbitrarily 
many spectral modes of equal quality. However, as the spectral 
mode index increases, the radiation resistance tends to zero, 
which would make practical design of antennas harder. If we 



only use spectral modes with big enough radiation resistance, 
then we recover the result in [2], i.e., the degrees of freedom 
is proportional to the surface area of the source. 

The following notations are used: (■)*, (•)* and (-y are vec- 
tor complex conjugate, matrix transpose, and matrix conjugate- 
transpose. jn\-) and hn(-) denote the spherical bessel func- 
tion of the first and the third kind, Y nm (-, •) is the spherical 
harmonics. The unit vectors on the spherical coordinate are 
f , 6 and <fi. I„ is the n x n identity matrix. 

II. Radiation Power over Modes 

A. Known results on singular value decomposition 

Assume some current source exists in V £ K 3 , and not 
in V . Then, the dyadic Green function G(r,r') relates the 
source current to the electric field on V as 



E(r) = iu>n [ G(r, r') • J(r')dr', r e V c 

Jv 



(1) 



where J(r') represents the current source at r', uj is the carrier 
frequency, and /i is the permeability of V c |8, p. 376]. As 
shown in (2), the Green function can be decomposed into 



oo n 



G(r,r')^fc£ J2 E^Vl) U «™'( r ) V ^( r ^ 



n— 1 m— — n , 



n(n + 1) 

where k is the wave number and 

IW(r) = V x rh£\kr)Y nm (6,4>), 
U nm2 (r) = iv x V x rh^{kr)Y nm (fi,4>), 



(2) 



V„ TO i(r) = V x rj n (kr)Y, 



■urn y ' > f ) 



V nm2 (r) = -V x V x rj n (kr)Y n 
k 



Now, assume V € M 3, is the interior of a sphere with radius 
Rv and the region of field measurement S is a surface of 
sphere with radius Rs > Rv, where both spheres are centered 
at the origin as shown in Fig. [T] Defining normalized eigen 




transmitter 



receiver 



Fig. I. Illustration on scattering environment. 



functions and eigenvalues as 

/ \ U nm / (r) A 1 . . 

u nm i{r) = . = U nmi (r), (3) 

/j s \V nml (r)\ 2 dn ^mi 



v„ m i(r) = 



L|V„ m; (r)| 2 dr 



is \U nm i(r)\ 2 dSly/ f v \V nm i(r)\ 2 dr 
n{n + 1) 



the Green function can be rewritten as 

oo 

G(r,r') =ifc^a p u p (r)v;(r'), (4) 

P =i 

where the mode index is defined as 

p = 2(n(n + 1) + m - 1) + I 

for n 6 N, — n < m < n, I — 1,2, so that one-to-one 
correspondence between (n, m, I) and p exists. The 
be decomposed into (2) 



source can 



(5) 



(6) 



j( r ) = J p v p( r )' 

P =i 

where 

J p = J(r)-v p (r)dr,peN 
Jv 

B. Radiation power analysis 

Assume modes p = 1, ...,M are used. Using Q and ([6j, 
([TJ can be rewritten as 

A/ A/ 

E(r) = -u>[xkJ2 VpJp'Vpi*) / v p(0 ' V( r ') dr ' 

M 

= -7]k 2 ^a p J p Up(r),r e S, 



(7) 



P =i 



by applying the orthogonality of v p (r) and = r/k, where 
rj = 1207T is the wave impedance. 

Define a M x 1 channel input X = [Ji, Jm]*> a M x M 
matrix £ = diag {ui, cta/} consisting of singular values, 
and a M x 1 vector *(r) = [ui(r), ...,u A /(r)]*. Then, 

E(r) = -ryfc 2 (^(r)) 4 £X,r e 5. (8) 




Fig. 2. Radiation resistance over modes when Ry = 2A. 



Using the definition above, the following theorem can be 
derived. 

Theorem 1. The amount of power leaving the surface S 
generated from the field given in ([8]) is 

rjk 2 



-X f S 2 T 2 X 



where T = diag {l/y/T\, 1/,/tm} and 



(9) 




h%(kR s ) 



where n is a corresponding index for each p = (n, m, I). 
The proof is in Appendix [A] 

Definition 1. The radiation resistance R r 
p = 1,2, is defined as 

r)k 2 <j 2 p 



for each mode 



2R v t p 



[Ohm] 



Note that the radiation resistance is irrelevant to R$- It 
decreases as n increases as shown in Fig. [2] 

Remark 1. The radiation power on mode p = 1,2,..., can be 
represented as P p = R ra d,p \Jp\ 2 Rv- 



III. System Model 

A. Power constraint 

Suppose signals are transmitted via different modes over 
h channel uses, and the source exists in V as shown 



in Fig. 3(a) Also, assume codewords are x"(m) = 
(xi(m), ... 7 Xn(fh)) , m € {l,...,2 nR } where R is the trans- 
mission rate, and x^(m),m € {l,...,2" R }, is the channel 



II-B 



input at time i £ {l,...,n} based on the definition in 
Applying Theorem [j] an average radiation power constraint 
can be defined as 



■qk 2 



^x](m)S 2 T 2 Xi(m) < nP, fn £ {l,...,2 flR } . (10) 




(a) 



(b) 



Fig. 3. Illustration of how receiving antennas work, (a) Receiving short 
dipole antennas uniformly distributed on S (b) Equivalent Thevenin circuit 
for short dipoles. The left half models the antenna, while the right half is the 
load. 



B. Receiver 

Based on the assumptions in Section [II] the electric field 
generated by the current in V carries the information traveling 
through S. In wireless communication, after dispersion and 
reflection, it finally arrives at the receiver as shown in Fig. [T] 
where the receiver is confined in a spherical region Vr and 
Sr is a sphere surrounding Vr. 

The overall channel has three parts, from V to S, from S to 
Sr, and from Sr to Vr. For simplicity, we assume the receive 
antennas (short dipole antennas) are located on S instead in 



Vr as shown in Fig. 3(a) and characterize the capacity of the 
channel from the source in V to the receive antennas on 

Let s q = (Rs, q , 4>q), q = 1, ■ ■ . ,N denote the location of 
the g-th dipole. We assume N 3> 1 for sufficient averaging 
and dipole locations are uniformly chosen on S. Assume the 
length of the dipole is L and e g denote the orientation of the 
q-th dipole for all q — l,...,N. 

Each short dipole can be modeled by an equivalent circuit as 
depicted in Fig |3(b)| where Vr is the voltage coming from the 
incident wave, R r is the radiation resistance of the dipole, R^ 
is the loss resistance, Rt is the resistance of the load, Xa is 
the antenna reactance, and Xt is the reactance of the load [9 
p. 84]. We choose R L = 0,X A = -X T and R r + R L = R T 
for conjugate matching to deliver the maximum power to the 
load. Then, the power transmitted to the load resistor is 



Pt = 



SRt 



C. Channel construction 



Define the receive signals without noise as 



V Tq 



,q=l,...,N, 



(11) 



so that its square is the received power by the load of the q-th 
antenna. In addition, 

2 



Rt„ — Rr„ — 80 



irL 



,q=l,...,N, 



(12) 



'Note that the path from Sr to Vr can be analyzed similarly as that from 
V to S and the path from S to Sr is a simple scattering environment. This 
gives the motivation for the simplification. 



Due to the uniform current density, we have 

V Tg = E q L,q = l,...,N. (13) 

where 

E q =E(s g ) -e q ,q = l,...,N, 

assuming the incident field is a plane wave |9] p.91]. Including 
the thermal noise across Rt, the received signal can be 
rewritten as 

V Tq 



Vq 



\/8RTq 

A 




E(s g ) • e 9 + z q 

OC 

2J (Up(s g ) ■ e g ) <T P J P + Zq ,q= 1, N, 
P =i 

which follows by using d7]l. Using the matrix representation, 
we have a Gaussian vector channel given as 




$sx + z 



where Y 



,2/jv]',* = [<t>x,...,<t> 



Ml >' 



[Zl,- 



(14) 

,z N \ 



such that Mp = 1,2, ...,0 p = [u p (si) • e 1: ...,u p (sjv) • ejv] ■ 
We assume Z ~ CM (0, Nq WTjy) where N is the noise 
spectral density, and W is the channel bandwidth. 

IV. Capacity Upper Bound 

Assume a total of N points are uniformly sampled on S, 
where N is assumed to be a multiple of 3. We assume the 
orientation of each dipole is given as 

r, q = 1,4,..., N -2, 
0, q = 2,5,...,N-l, 
4>, q = 3,6, N. 

A. Approximated channel capacity 

Using the channel given in dl4}, the channel capacity is 

1 



C = max W log 

Kj:tr(K s )<P 

= max W log 

K s :tr(K i )<P 



= max W log 

K jt :tr(K i )<P 



N w 
1 3tt 



N w 2 
1 3tt 



GK x G f 

*T - 1 K*T _1 * t 



T"i$t$ x i K 



U' 



(15) 



(16) 



where p5| ) follows since = ^TSKxST by assuming 
TEX,N W = N W, and ((T6j) follows since 




M 



BAI. Also, we have 



N 



9=1 



Using <|3j and basic calculus, we have 
u*(r) -u p/ (r)dn 



N 



Jim y){%(s 9 )-e 9 }*{v(s,)-ej 

AT— foo * — ' 



9=1 



47T 

(JV/3) 



s PP ',yp,p' = 1,2, 



Since 



4>l<t>p> = j^Spp' +o(N),\Jp,p' = 1,...,M, 
continuing from ( |16) , we get for sufficiently large AT, 



C = max log 

K r MK % )<P 



1 iV\ 



iV, 



-T Ki 



w 



(17) 



(18) 



B. Upper bound 

Since the field is measured based on the plane wave 
assumption, M should be appropriately chosen so that each 
dipole antenna receives an asymptotic plane wave, even if 
the wave itself is spherical. Note that this restriction on M 
does not hurt the capacity since we later let R$ tend to 
infinity and thus M can also tend to infinity. In order to 
satisfy the condition, the incident wave should asymptotically 
be transverse electromagnetic (TEM) wave near the receiving 
region. In |jT0} p.925], 

p ikRs ™ 



(n + m)! 



= H) 



n+l ' 



kR s 2 m m\(kR s ) m (n - m)\ 

a ikRs 



kR s 



(19) 



where ( fT9| i follows when kRs ^> n 2 , and it is sufficient to 
satisfy the TEM condition. 

Definition 2. Assume the total radiation power is P and the 
total power received on the antenna is Pr. Then, total effective 
sold angle fl e is 

Pr 



Vl P = 47T X 



(20) 



Theorem 2. With a total radiation power constraint P, the 
capacity is given as 

aP 

C= N lOg€ > 
where a = 5h-,0 < a < 1. 

Proof: The sum of the power received on dipoles is 

3N 2 , 2 M M N 

Tj rC 



we can choose 



N = [8afc 2 i?|] , 



(23) 



in order to separate the dipoles to prevent the coupling effects. 
For a given sufficiently large Rs, we can choose Uq such 



that n 2 , < kR s . Let M = (n - l)(n + 1), which is the 
number of modes satisfying n < no — 1. Then, we have 



1 



k 2 R 2 



1 + 



(24) 



\kR S/ 

Now, assume 8ak 2 R s is an integer for simplicity. Let tiq = 
koy/Rs when fco "C k, and then, M — k 2 Rs — 1. By using 
(1531 and dM , 



C = max W log 

K i :tr(K i )<P 



max W log 

%:tr(K*)<P 



Im + 



= MWlog ( 1 
aP, 



Im 

MJVw 



N w 
a K- 



(25) 



as M — > oo provided that i?s — > oo. Note that we only require 
Rs — > oo due to a limitation of our analysis (plane wave 
assumption for short dipole antennas). We believe our result 
would hold for any finite Rs > Rv- ■ 

Appendix A 
Proof of ^ 

For notational simplicity, if inputs are omitted, 
hn\Y nm ,P nm indicate h£' (kr), Y nm (6,<t>), P„ m (cos 6), for 
all n and m. Both electric and magnetic field on S satisfy 
the Maxwell equations on the source-free region such that 

V x E = iw/zH, (26) 
VxH = -icjeE, (27) 



where e is the permittivity on S (8J. Using ( |26| l and ( |27| i, we 
have 



V x V x E = fc 2 E, 
From Q and po*] ), the magnetic field is given as 

H(r) 



(28) 



ik / <7pJp V x u p (r). 



/ , u p""p 
p=i 



g=l 



Using (40) and (42) in (2), we know 

U„ m i(r) = VY nm x rrh£\ 

i 



160 



J2 u *p( s i) •vKX ct p' J p j p' U„ m2 (r) = fn(n + 1)^-F r 



o=l p'— 1 q=l 



kr 



dh 



dr 



?7 2 fc 2 Ar 
1920tt 



A/ 



2 

Pi ) 



(21) 



p=i 



Also, (41) and (44) in (2) show 

= n(n + l) h^(kR s ) 



for sufficiently large AT, where pTj ) follows by ( fTTj l. Since 
both transmit and received power are the same in case 



(1* 

r 2 



AT = \8k 2 R 2 s ] 



(22) 



n(n + 1)* 
2n+ 1 



4 1 2 1 (fci?5) 



n 2 (n + 1) 
2n+ 1 



^+i(^s) 



The complex power flow leaving the surface S is 

P c = <f ^(E(r)xH*(r)).ds 

j s 



ipp' 



(29) 



p— 1 p'==i 

where 

7 PP ' — j> Up (r) x V x u*,(r) • ds 

1) p = (n, to, l),p' = (n'mf, 1): 
Cp Cp' 7pp' 

= fc^ U nml (r) x t£, m , a • rr 2 dfi 
= fc J U nm i(r) • U*, ro , 2 x rr 2 dfi 



(30) 

VF„ m x r • VY*„„ x rdfi 



Cp Rs + R% 



i dfclV 



hi) 



(!) dr 



&nn'&mm'i (31) 



r=Rs 



where ( |30| > follows by vector identity (a x b) • c = a • (b x c), 
and pT| follows since 

n(n + 1) 



VY nm x r • VY* W x rdn 



which follows by (9.120) in (TT). Thus, 

2) p — (n, m, = (n'mf , 2); 



7pp> oc y U„ mi (r) x U*, ma (r) ■ ds 

cx y (VY nm x r) x (VF„* m , x f) • ds 

« J Vy n* m' • (r X Vr «m) ««1 

cx y f x VF n *, m , • r x (f x VY nm ) dft 
= 0, 



(32) 



(33) 



where ( |32| ) follows by vector identity a x (b x c) = (a • c)b — 
(a • b)c, and ([33} follows by (9.121) in jn]. 

3) p — (n, to, 2),p' = (n'mf , 1): 

7pp, cx y U„ m2 (r) x U;, m , 2 (r) ■ fdQ 

cx y r x U nm2 (r) • U;, m , 2 (r)dn (34) 
= 0, (35) 



where the proof is done similarly as in A-2 



4) p = (n, to, 2),p' = (n'mf , 2): 
CpCp'Tpp' 

= kjrx U nm2 (r) • U,;, m , 1 (r)r 2 ^ 



<//) 



(i) 



r/r 



vr„ m x f • vr, m , x rdn 



1 dh u 



r=R s ■ 



where the proof is similar with ( f30| > and ( (3T| . Thus, we get 



7pp' 



r 2 



C 2 



l d/il ir 



(36) 



r=Rs 



Note that power of each mode is independently generated 
since 7 pp / = when p ^ p'. Assuming 



1. 



p = 1,3, 
p = 2,4,..., 



the radiation power is given as 
5He{P c } 

= ^g Tpa2|Jp|23m . 

p— i ^ 

where ( |37| i follows by using 



1 dh { n ] 



r=R s ■ 



(37) 



Jm ■ 



1 ctt, 



1 



k 2 R s \h£\kRsW 



r=Rs 
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